A Bright Spatially- Coherent Compact X-ray Synchrotron Source 
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Each successive generation of x-ray machines has opened up new frontiers in science, such as 
the first radiographs and the determination of the structure of DNA. State-of-the-art x-ray sources 
can now produce coherent high brightness > keV x-rays and promise a new revolution in imaging 
complex systems on nanometre and femtosecond scales. Despite the demand, only a few dedicated 
synchrotron facilities exist worldwide, partially due the size and cost of conventional (accelerator) 
technology [1 . Here we demonstrate the use of a recently developed compact laser-plasma acceler- 
ator [21 El [4] to produce a well-collimated, spatially-coherent, intrinsically ultrafast source of hard 
x-rays. This method reduces the size of the synchrotron source from the tens of metres to centimetre 
scale, accelerating and wiggling a high electron charge simultaneously. This leads to a narrow-energy 
spread electron beam 02 2] an d x-ray source that is > 1000 x brighter than previously reported 
plasma wiggler [5] |6] and thus has the potential to facilitate a myriad of uses across the whole 
spectrum of light-source applications. 



I. INTRODUCTION 

There are many proposals to use extreme non-linear 
interactions of the latest generation of high-power ultra- 
short pulse laser systems to produce beams of high energy 
photons with high-brightness and short pulse duration. 
For example, high-order harmonic generation promises 
trains of coherent pulselets [3 |8] whilst Compton scat- 
tering could extend energies into the 7-regime [9l HO] . 
However, these and other similar schemes are complex 
and have stringent laser requirements. An alternative 
proposal has been the use of compact laser-plasma accel- 
erators to drive sources of undulating/ wiggling radiation. 

These accelerators use the plasma wakefield generated 
by the passage of an intense laser pulse through an un- 
derdense plasma [11]. Such wakefields can have intrinsic 
fields > 1000 x greater than the best achievable by con- 
ventional accelerator technology, and thus can accelerate 
particles to high energies in a fraction of the distance. Re- 
cently, it has been demonstrated that at high laser power, 
the wakefield can be driven to sufficient amplitude that 
it can trap particles from the background plasma and ac- 
celerate them in a narrow energy spread beam[2j [3j [4], 
now producing beams of electrons of ~ GeV energy on 
the order of 1 cm (H [13] . 

Such electron sources are clearly of interest to replace 
the accelerators which drive current synchrotron sources, 
and typically use multiple periods of alternately poled 
magnets (undulators or wigglers) to reinforce the syn- 
chrotron emission over a length of ~ metres. The first 
demonstrations of wakefield driven radiation using exter- 
nal wigglers have now indeed been reported, though still 
being limited to optical or near-optical wavelengths and 
modest photon numbers [T4j [15] . 

However, the particles being accelerated in the plasma 
accelerator also undergo transverse (betatron) oscilla- 



tions due to the focusing fields of the plasma wave. 
The oscillations occur at the betatron frequency oup = 
cjp/y^, where co p is the plasma frequency and 7 is the 
Lorentz factor of the electron beam. Since this beta- 
tron wavelength (A^ « 27rc/ujp) is much smaller than 
the period of comparable external wigglers, again due 
to the extremely large electric fields of the wakefield, 
the wavelength at which these particles radiate (of order 
^ ~ A/3/7 2 ) can be in the x-ray region for 7 > 200 i.e. for 
an electron beam of energy E > 100 MeV. For small 
transverse amplitude and thus small strength parameter 
K = jrpcup/c <C 1 (i.e. undulator limit), the radiation 
spectrum will be narrowly peaked about a fundamental 
energy. As K — > 1, radiation also appears at harmon- 
ics. For large K ^> 1 (large amplitude - wiggler limit) a 
synchrotron spectrum with broad emission consisting of 
closely spaced harmonics is produced [TO] . 
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Here, / is the radiated energy, Np is the number of oscil- 
lations and E CT - lt = ZhK^uop is twice the energy above 
and below which roughly half of the total energy is ra- 
diated. /C2/3 is a modified Bessel function of the second 
kind. For E > E cr it, the spectrum of radiated energy 
decays exponentially. The radiation is confined to a cone 
with opening angle ~ K/j z - 

Previous measurements of betatron radiation from 
a broad energy spread electron beam from a laser 
wakefield accelerator have demonstrated x-ray emission 
with energies up to ~ 1 keV with a brightness up to 
10 19 ph/s/mrad 2 /mm 2 /0.1%BW [5 . Here we demon- 
strate the use of a compact accelerator (in total < 1 cm), 
which produces a beam with much higher charge at high 
energy (~ 250 MeV), to produce a 1-100 keV x-ray source 
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with more than 1000 x greater brightness than previous 
laser driven betatron sources [5j [6]. Furthermore, we 
demonstrate that this radiation source exhibits evidence 
for spatial coherence. 



II. EXPERIMENTAL RESULTS 

The experiment was performed by focusing an intense 
short pulse (~ 30 fs, ~ 2 J) laser onto the front edge of 
3, 5 and 10 mm helium gas jets (see methods). Electrons 
beams with narrow energy spread features were observed 
from all nozzles, at densities of 4 — 22 x 10 18 cm -3 . Due 
to consecutive phases of injection [13], the electron beam 
consists of multiple filaments. For example, for a density 
of 8 x 10 18 cm -3 on the 5 mm nozzle, electron beams 
of E = (230 ± 70) MeV with AE = (25 ± 10)% energy 
spread at FWHM were observed with an average of 2.2 ± 
0.4 filaments per shot, with a root-mean-square (RMS) 
divergence of 1.5 x 1.8 mrad 2 . The beam position had a 
RMS pointing fluctuation of 4.8 x 4.7 mrad 2 . 

The average and maximum energy of the electron beam 
follow the typical wakefield density scaling law [I3j [17] . 

With the electron beam deflected away from laser axis 
by the spectrometer magnet, a bright (undeviated) beam 
of x-rays was also observed co-propagating along the laser 
axis. It was imperative to first prove that this x-ray 
source originates from the plasma itself. In order to do 
this a grid of 60 /im diameter silver wires was placed a 
few centimetres from the target. X-rays originating from 
the interaction region project the outline of the mesh 
onto an imaging plate [18 . A strongly collimated beam 
of x-rays is evident in figure la. When either the laser 
power or plasma density was reduced to inhibit the elec- 
tron beam, the x-ray beam also disappeared, showing 
that the generation of the x-rays is linked to the electron 
beam. The profile is elliptical with a FWHM divergence 
of X x y = 4 x 13 mrad 2 , corresponding to a wiggler 
parameter K = #7 of K x = 1.5 and K y = 5 for a simul- 
taneously measured electron beam energy E = 200 MeV. 
The x-ray beam pointing is extremely stable, as can be 
deduced from figure lb, which shows the sum of five con- 
secutive shots. Their combined divergence is not signifi- 
cantly larger than that of a single shot measurement, the 
RMS pointing stability is 5 mrad in the horizontal and 
vertical direction, similar to the pointing stability of the 
electron beam. 

To give an indication of the x-ray source size, micro- 
scopic objects were backlit with the x-ray beam. Figures 
lc,g show photographs of collections of wires of various 
sizes, and figures 1 d,f,g,h show the corresponding radio- 
graphs of these objects. Even features as small as a 5 
/am wire are resolved, indicating that the betatron x-ray 
source is < 5 /im, smaller than the size of the plasma wave 
in which the radiating electrons were trapped (« 20 /im 
diameter). 

To quantify the source size more precisely, a half-plain 
was backlit with the x-ray beam. A typical intensity dis- 
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Figure 1 



FIG. 1: X-ray beam profile behind a wire mesh from a single 
shot (a) and five consecutive shots (b) showing a 5 x 11 mrad 2 
beam with 5 mrad pointing stability. Photograph (c) and x- 
ray radiograph (d) of a 12 /im thick gold foil with rectangular 
500 to 20 /xm features. Photograph (e) and x-ray radiograph 
of wire triplets (f: 50/250/100 /im), (e,g: 50/100/20 /im), (h: 
5/20/10 /im). All features are resolved indicating a betatron 
source < 5 /im. (i) Single edge Fresnel diffraction from a 
cleaved InSb crystal produces a half-shadow, consistent with 
a (2.0 dz 0.5) /im spatially coherent source. 



tribution on the detector looks like a half-shadow (figure 
li), whose details convolve information about the x-ray 
source and half-plain. The half-plain was a 0.25 mm thick 
cleaved InSb crystal (< 6% transmission below 20 keV) 
and resembles an ideal step function. Assuming a gaus- 
sian source function, convolution with the step function 
gives an error function which fits the experimental in- 
tensity distribution well. Figure 2 (inset) shows several 
experimental intensity lineouts and fitted error functions 
which predict the sharp rise and provide an estimate of 
the source size, in the range 2.9 - 6.0 /im at FWHM for 
the traces shown in figure 2 (inset). The fit with the er- 
ror function, however, fails to reproduce the ringing that 
follows the sharp rise on the high intensity side. 

To accurately model the shape of the intensity dis- 
tribution, it is necessary to use Fresnel diffraction (see 
methods) , where the details of the diffraction pattern de- 
pend on the spatial and spectral intensity distribution of 
the source and the dimensions of the setup. Figure 2 
shows an experimental and modelled intensity distribu- 
tions based on a gaussian spatial intensity profile and 
synchrotron spectrum, assuming a constant phase across 
the source. For the solid red curve in figure 2, E cr i t = 8 
keV and an 1/e 2 intensity radius of w x = 1 /im was as- 
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sumed, which best reproduced both the sharp rise and 
the amplitude and width of the first fringe. Changing w x 
or E crit under or overestimates the height and/or width 
of the overshoot and/or the rise, as shown by the other 
curves in figure 2. This determines that the 1/e 2 radius of 
the source is in the range 0.5 and 2 fim and E cr a between 
4 and 16 keV. 

The ringing must be due to interference of the radia- 
tion that has originated from different regions of a spa- 
tially extended source. If the phase relationship across 
the source was random, the fringes would be blurred, 
as shown in gray open circles in figure 2. The spectral 
width of the source also reduces fringe visibility. The 
fact that the data exhibits a clear first fringe can only be 
attributed to spatial coherence of the source. 

The x-ray yield measured with an x-ray CCD was stud- 
ied as a function of electron charge. The x-ray signal 
increases with interaction length from 5 mm to 10 mm 
(Figure 3a). This is not surprising since a longer in- 
teraction length can accommodate more betatron oscil- 
lations. However, the x-ray yield per charge increases 
more than linearly from the 5 mm to 10 mm nozzle, by 
an average factor of 3. The electron beam exhibits os- 
cillatory transverse deflections in the spectrometer nor- 
mal to the dispersion direction, which is evidence for the 
presence of betatron oscillations during acceleration [19] . 
Indeed a high x-ray yield is correlated with the promi- 
nence of these betatron remnants (Figure 3a inset). An 
interaction of the accelerating electron beam with the 
laser (betatron resonance), can promote an increase in 
oscillation amplitude [6 . For this to occur, the elec- 
tron beam must catch up with the laser, which would 
be favoured by a longer interaction. The faster than lin- 
ear scaling of the x-ray yield may be a coherent effect 
based on micro-bunching, which simulations have shown 
as a consequence of a betatron resonance [20] . 

The spectral properties of the betatron radiation were 
determined by measuring the x-ray transmission through 
a set of filters with an x-ray CCD detector [6 . Assuming 
the spectrum is synchrotron-like, the measured E CY - lt in- 
creases with density and electron energy. For the 5 mm 
nozzle at n e = (1.0 ± 0.4) x 10 19 cm -3 , the measured 
^crit = (29 ± 13) keV as shown in figure 3b. 



III. NUMERICAL SIMULATIONS 

Numerical modelling was carried out using electron 
trajectories obtained from the fully relativistic PIC code 
OSIRIS. The simulations were run in the boosted frame, 
which due to relativistic length contraction and time 
dilation, offer higher resolutions for shorter run times. 
The trajectories were post-processed to yield the char- 
acteristics of the betatron radiation (see methods). Fig- 
ure 4a shows a simulation for a 5 mm nozzle at n e = 
8 x 10 18 cm -3 , the spectrally integrated x-ray beam pro- 
file with a divergence of 5 x 16 mrad 2 . The profile is 
elongated in the direction of laser polarisation, demon- 




-0.2 -0.1 0.1 0.2 0.3 0.4 
x/mm on camera 

Figure 2 

FIG. 2: The x-ray source casts a shadow of a half-plane on 
the detector, (inset) Intensity distribution of the half-shadow 
integrated along the edge (solid black) and fit with error func- 
tions (solid grey) for four different shots, giving source sizes 
of 2.9, 3.8, 6.0 and 3.9 fim at FWHM from top to bottom, 
assuming gaussian source distribution. Error functions mimic 
the observed intensity distribution closely but fail to repro- 
duce the oscillations, (main) Close-up of measured (open 
black squares) and modelled intensity distribution using Fres- 
nel diffraction methods, assuming a spatially coherent gaus- 
sian source with synchrotron spectrum E CY i t /wo of 8 keV/1 
/im (solid red), 8 keV/3 /mi (dashed green) and 4 keV/1 /xm 
(dash-dotted blue) and an incoherent source with 8 keV/2 /im 
(open grey cicles). 



strating excellent quantitative agreement with the mea- 
surement (figure la). Figure 4b shows the modelled x-ray 
spectrum as a function of energy and angle. The x-ray 
flux peaks on-axis at 10 keV and extends to ~ 100 keV. 
A lineout taken on-axis, as would be measured by our 
detector, is plotted in figure 3b, comparing well with 
the measured spectrum. A total of 10 8 photons are pre- 
dicted between 1 and 84 keV from a 30 fs electron bunch, 
where we measure 10 6 to 10 8 , depending on electron 
charge, corresponding to a peak brightness 10 21 — 10 23 
ph/s/mrad 2 /mm 2 /0.1%BW. The x-ray pulse duration is 
of the same order as the electron bunch duration, and 
is thus of ~fs timescale. The modelling reveals that the 
fastest electrons oscillate with the smallest amplitudes 
rp < 2/im, consistent with the measured source sizes. K- 
parameters range from 1-30 with an average of 7 close 
to the measured experimental value. 



IV. CONCLUSIONS 

In conclusion, we have shown, that a laser driven 
plasma can simultaneously serve as particle accel- 
erator and wiggler, producing high quality beams 
of x-rays, with mrad divergence, fim source size, 
tens of keV critical energy and peak brightness ~ 
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FIG. 3: (a) X-ray yield as a function of electron charge, for 
the 5 mm (black squares) and 10 mm (grey triangles) noz- 
zle. The x-ray signal increases faster than linearly with in- 
creasing interaction length (and thus Np). (a inset) shows 
an energy-dispersed electron spectrum with prominent trans- 
verse features remnants of betatron oscillations, (b) The be- 
tatron spectrum obtained from numerical modelling agrees 
well with the experimentally measured spectrum within one 
standard deviation. 



10 23 ph/s/mrad 2 /mm 2 /0.1%BW. The measured radia- 
tion is bright, spatially-coherent, and intrinisically ul- 
trafast opening up a multitude of advanced applica- 
tions, such as phase contrast and lensless imaging, previ- 
ously only possible with large conventional light sources 
[2T| [22] . The laser plasma wiggler has the potential 
of making high-brightness radiation sources widespread, 
and thus impacting all areas of science and technology. 
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Figure 4 

FIG. 4: Numerical modelling of interaction; (a) spectrally in- 
tegrated x-ray beam profile showing an elliptically elongated 
beam profile in the direction of laser polarisation, (b) angu- 
larly and spectrally resolved x-ray flux displaying a peak on 
axis at 10 keV and a tail extending to ~ 100 keV. 



with a central wavelength of Ao = 800 nm, a gaus- 
sian full-width half maximum (FWHM) pulse duration 
of r = 32 fs and a maximum energy of 2.3 J were focused 
to a peak intensity of I = 4.7 x 10 19 Wcm -2 or normalised 
vector potential ao = 4.7 with an off- axis parabolic mir- 
ror of focal length / = lm and F-number of F = 10. The 
transverse intensity profile of the laser in vacuum yields 
a focal spot diameter of oVwhm = 11 -2 /im. Typically 
55 % of the pulse energy is within oVwhm- 



Presnel Diffraction 



METHODS 



A. Laser 



The experiments were carried out on the Hercules laser 
at the University of Michigan. Linearly polarised pulses 



Fresnel diffraction occurs when a wave passes through 
an aperture and diffracts in the near field, causing any 
diffraction pattern observed to differ in size and shape, 
depending on the distance between the aperture and the 
projection z, the size of the aperture A and the wave- 
length of the wave A. In the present case, i 1500 mm, 
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A < 5 x 10 -7 mm (> 2 keV) and a ~ 125 fim assum- 
ing an x-ray beam of 5 mrad divergence and a half-plain 
aperture at o = 50 mm from the source. Therefore the 
Fresnel number F = A 2 / (iX) « 16 > 1, requiring a treat- 
ment in the Fresnel regime. The intensity distribution of 
a monochromatic point source on the projection is given 

by 

/(*) = E-E* = | | Q + C(w(x))j 2 + Q + S(w(x))y 

(2) 

where E and E* are the electric field and its complex 
conjugate, C and S are the Fresnel functions, w(x) = 

x (o/(i + 6) ^/2(l/z + l/o)/X^j and x is the position on 

the projection. 

If the source is not point-like, the electric field has to 
be convolved with the source function and spectral dis- 
tribution E(\, x) = J J x R(X)g(y)E(X, x — i/o • y) dydX 
to compute I(x), where R(X) is the spectral response of 
the detector. 



which a linearly polarised pulse with experimental pa- 
rameters is focused 0.25 mm into the plasma. The 
plasma density profile increases linearly from zero to 
n e = 8 x 10 18 cm -3 in the first 0.5 mm, it is constant 
for 3.5 mm, and falls linearly to zero in 0.5 mm. Simula- 
tions were performed in a relativistically boosted frame 
(7 = 5), which allows for significant computational gains 
[24] . The simulation box corresponds to 60 x 94 x 94 /im 3 
in the laboratory frame. A total of 9.2 x 10 8 particles 
were pushed for 5 x 10 3 iterations. The resolution in 
the laser propagation direction z is koAz = 0.11, and 
kpAx = kpAy = 0.16 in the transverse directions. To ob- 
tain the radiation emission from the simulation, a post- 
processing diagnostic is used on a set of trajectories of 
injected electrons [25]. The position and momentum his- 
tory of the electrons is used to deposit the radiated fields 
on a virtual detector. 



C. Numerical Modelling 

Three-dimensional numerical simulations were per- 
formed with the particle in cell code OSIRIS [23], in 
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